home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / CHPDI.z / CHPDI
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. CCCCHHHHPPPPDDDDIIII((((3333FFFF))))                                                            CCCCHHHHPPPPDDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CHPDI   - CHPDI computes the determinant, inertia and inverse of a
  10.      complex Hermitian matrix using the factors from CHPFA, where the matrix
  11.      is stored in packed form.
  12.  
  13.  
  14. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  15.       SUBROUTINE CHPDI(AP,N,KPVT,DET,INERT,WORK,JOB)
  16.  
  17. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  18.      On Entry
  19.  
  20.      AAAAPPPP COMPLEX (N*(N+1)/2)
  21.         the output from CHPFA.
  22.  
  23.      NNNN INTEGER
  24.         the order of the matrix A.  KVPT    INTEGER(N)
  25.         the pivot vector from CHPFA.
  26.  
  27.      WWWWOOOORRRRKKKK COMPLEX(N)
  28.         work vector.  Contents ignored.
  29.  
  30.      JJJJOOOOBBBB INTEGER
  31.         JOB has the decimal expansion  ABC  where
  32.         if  C .NE. 0, the inverse is computed,
  33.         if  B .NE. 0, the determinant is computed,
  34.         if  A .NE. 0, the inertia is computed.
  35.         For example, JOB = 111  gives all three.  On Return Variables not
  36.      requested by JOB are not used.
  37.  
  38.      AAAAPPPP contains the upper triangle of the inverse of
  39.         the original matrix, stored in packed form.
  40.         The columns of the upper triangle are stored
  41.         sequentially in a one-dimensional array.
  42.  
  43.      DDDDEEEETTTT REAL(2)
  44.         determinant of original matrix.
  45.         Determinant = DET(1) * 10.0**DET(2)
  46.         with 1.0 .LE. ABS(DET(1)) .LT. 10.0
  47.         or DET(1) = 0.0.
  48.  
  49.      IIIINNNNEEEERRRRTTTT INTEGER(3)
  50.         the inertia of the original matrix.
  51.         INERT(1)  =  number of positive eigenvalues.
  52.         INERT(2)  =  number of negative eigenvalues.
  53.         INERT(3)  =  number of zero eigenvalues.  Error Condition
  54.  
  55.      AAAA division by zero will occur if the inverse is requested and  CHPCO  has
  56.      set RCOND .EQ. 0.0 or  CHPFA  has set  INFO .NE. 0 .  LINPACK.  This
  57.      version dated 08/14/78 .  James Bunch, Univ. Calif. San Diego, Argonne
  58.      Nat. Lab.  Subroutines and Functions BLAS CAXPY,CCOPY,CDOTC,CSWAP Fortran
  59.      ABS,CABS,CMPLX,CONJG,IABS,MOD,REAL
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.